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Abstract 

Since the 1990s chaotic cat maps are widely used in data encryption, for their very complicated 
dynamics within a simple model and desired characteristics related to requirements of cryptography. The 
number of cat map parameters and the map period length after discretization are two major concerns 
in many applications for security reasons. In this paper, we propose a new family of 36 distinctive 3D 
cat maps with different spatial configurations taking existing 3D cat maps [l]-[4] as special cases. Our 
analysis and comparisons show that this new 3D cat maps family has more independent map parameters 
and much longer averaged period lengths than existing 3D cat maps. The presented cat map family can 
be extended to higher dimensional cases. 

Index Terms 

Arnold Transform, Image Encryption, Cat Map, Automorphism 

I. Introduction 

In the last decade, many efforts have been recognized to study behaviors of dynamic systems and related 
applications. As one type of dynamic systems with very complicated behaviors, chaos systems are widely 
reported in mathematics, physics, engineering, economics, and biology. Especially for cryptography 
and encryption, chaotic cryptosystems [5] are demonstrated to have many analogous to conventional 
encryption methods [2], e.g. a conventional algorithm is sensitive to keys, while a chaotic system is 
sensitive to its initial values and parameters. Chaotic encryption systems also provide better solutions 
for image encryption where the conventional method are not suitable due to many intrinsic features of 
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images [2], e.g. bulk data nature and high correlations between pixels. Consequently, numerous chaotic 
encryption algorithms have been proposed in literature based on various principles [2], [6]-[8]. Among 
these applications, Arnold's cat map and its variants are commonly employed as a fundamental building 
block for data encryption [l]-[4], [7], [9]-[ll], watermarking [12], and pattern recognition [13]. 

Arnold's cat map [2], [14] is a discrete chaotic map named after Vlamimir Arnold. Specifically, it is 
erogdic and mixing, a C-system, a K-system and a Bernoulli system [14]. Mathematically, Arnold's cat 
map is defined in Eq. (1) 
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where the transform matrix C is defined as 
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In order to achieve higher map randomness and thus a more secure cryptosystem, the cat matrix C 
of Arnold's cat map is commonly replaced by Eq. (3), where new parameters help increase key spaces 
to resist brute-force attacks and longer averaged period lengths help improve map randomness to resist 
statistical attacks. 
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Based on Eq. (3), many efforts [l]-[4] are put on designing 3D cat maps to further increase the parameter 
space and map randomness. However, we notice that 1) they failed to consider all possible spatial 
configurations of a 3D cat map; and 2) their 3D cat map matrix elements are designed in more correlate 
rather than independent to each other. 

In this paper, we propose a new family of 3D cat maps with improvements on existing 3D cat maps 
[l]-[4]. This new family contains more parameters with longer averaged period lengths and taking these 
commonly used 3D cat maps as special cases. The rest of the paper is organized as follows: Section II 
reviews four existing 3D cat maps; Section III introduces our algorithm for generating the new 3D cat 
map family; Section IV analyzes the averaged period lengths of the proposed 3D cat maps with other 
3D cat maps; and Section V concludes the letter. 
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II. Existing 3D Cat Maps 

Based on the parametric 2D cat map in Eq. (3), Lian et al. [1] proposed a parametric 3D cat map 
C| D in 2003 by extending 2D cat map in zx plane and yz plane. 
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Later on, in 2004 Chen et al. [2] proposed a family of parametric 3D cat maps defined in Eq. (5) by 
composing three fundamental parametric 3D cat maps on xy and yz and zx planes defined in Eqs. (6), 
(7) and (8), respectively. 
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Based on Chen's cat map, Liu et al. [3] proposed an improved 3D cat map Cff in 2008 by introducing 
new parameters c and d as shown in Eq. (9) 
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Recently in [4], Pan et al. introduced their 3D cat map Cp D based on Chen's method, but in a different 
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formulation from Liu et al. 's. 
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It is worthwhile note that all these cat maps satisfy the condition that the determinant of 3D cat matrix 
is 1, i.e. 
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III. A New Family of 3D Cat Maps 
One common feature of existing 3D cat maps [l]-[4] that they are all extended from the 2D cat map 
Cpara defined in Eq. (3), which is indeed a general form of 2D cat but only one out of four possible 
spatial configurations as shown in Eqs. (11)-(14): 
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As can be seen, these configurations are inequivalent to each other in general because none of them has 
the unitary element located in the same position. However, spatial configurations are fail to be considered 
in previous 3D cat maps [l]-[4]. 



Meanwhile, we also notice that constructing 3D cat maps by directly extending the 2D cat map C 
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as done in [l]-[4] is not necessary, because all what we need is to construct a 3 x 3 matrix C 3D with 
the constraint that its determinant is 1 [2]. In other words, we are looking for a 3 x 3 matrix with 
the symbol set 6 = {a, b, c, d, e, f, g, h, i} as shown in Eq. (15) whose determinant satisfies Eq. (16). 
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We describe a general solution to the above constraint problem in Algorithm 1. This algorithm is able 
to generate 36 spatial configurations specified by parameter u and m for 3D cat maps, each configuration 
with six other independent parameters controlling cat matrix elements. 



Algorithm 1. Cat Matrix Generator 

Require: u is an integer in {1, 2, • • • , 9} 

Require: m is an integer in {1, 2, 3, 4} 

Ensure: C^? m is a 3D Cat matrix with |C^?J = 1 

1. Set symbol 6„ to 1 and compute the target determinant aei + bfg + cdb, - ug - bdi — aft) = 1 by substituting 

6„ = 1 

2. Form the four symbols in the positive or negative diagonal containing the unitary symbol & u to set $ 

3. Collect the terms containing the symbol $ m on the left side of the equation and leave all the other terms on 
the right side with a name rightside 

4. Set the coefficient term of symbol 3m to 1 and result in symbol $ m = rightside. 

5. Output C^? as a 3D cat map. 



Details of why this algorithm works are illustrated by the following example. Assume u = m = 1, 
then we first set the unitary symbol & u = &i = a = 1, i.e. we have 
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Consequently, the target determinant equation 




= 1 become 


ei+ bfg + corj 
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- bdi 
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We then obtain a symbol set $ = {e, i, f, h} for those either along the positive or the negative diagonal 
containing the symbol &i, and collect all terms containing symbol 3m = $i = t on the left side of the 
equation, leaving all the other terms on the right side, namely 

(i - cg)e = 1 - bfg - cOfj + bOi + ft) 

Next we set the coefficient term of symbol e to 1, i.e. i = cg + 1, and the simplified equation implies that 

e = l-bfg-cc)h + bc)(cg + l)+fh / 
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And we finish the 3D cat map construction, because in Eq. (17) is already of a 3 x 3 matrix with 

determinant 1. 

f l b c 

o to + ff> - bfg - cOf) + bcog + 1 f (17) 
f) eg + 1 

The other three variants with unitary element a = 1 of the 3D cat map family are shown below: 
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co + ei - ceg - boi + bcog - 1 
g bg + 1 i 

In general, a symbolic 3D cat map of the possible 36 spatial configurations can be easily obtained in a 
similar manner by feeding different us and ms in Algorithm 1, each configuration contains six independent 
parameters like those in Eqs. (17)-(20). Therefore, each new parametric 3D cat map is associated with 
eight independent parameters, two controlling spatial configurations and six controlling matrix elements. 
Consequently, this new family of 3D cat maps have more parameters but less correlated matrix elements 
than existing 3D cat maps [l]-[4]. Detailed comparisons about parameters and matrix elements between 
different 3D cat maps are summarized in Table I. 



TABLE I: Elements and Parameters in 3D Cat Maps 



Existing 3D Cat Maps 



Number of Items 












Constant Elements 


1 





2 








Constant 1 Elements 


2 





2 


1 


1 


1 Parameter Elements 


2 


2 


4 


4 


6 


2+ Parameter Elements 


4 


7 


1 


4 


2 


Parameters 


4 


6 


4 


4 


8 


Spatial Configurations 


1 


1 


1 


1 


36 
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It is worthwhile to note that the proposed 3D cat map family includes previous 3D cat maps [l]-[4] 
as special cases as shown in Table II, where each existing 3D cat map C 3D can be denoted by C^? 
with the eight parameters listed the table and symbol * indicates either constant elements or dependent 
elements on two or more parameters. For example, Cp 3 is a special case of Cffl um because Eq. (21) 
holds. 



i3D 



f~i3D |U=l,m=2 



b=a,c=c,t)=b,f=bc,g=d,l)=abcd 



(21) 



For verification, simply substitute these parameter values in C^? 2 . And we obtain * elements o = 1, 
e = bd + 1 = ab+1, and i = ccj+fh — bfg — coh + bcOg + l = cd+ (be) (abed) — abed— cb(abcd)+acbd+l = 
cd + 1, which are indeed the three corresponding elements of Cp D defined in Eq. (10). 



TABLE II: Denoting existing 3D cat maps with 
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Meanwhile, the new proposed 3D cat maps can be extended to a more general case simply by 
multiplying 3D cat maps with interested configurations. In general, we can construct such a mixed 3D 
cat map C^? as shown in Eq. (22), where s um is the indicator function for configuration (u, m) defined 
in Eq. (23). In this way, we are able to construct a more general 3D cat map with 36 new indicator 
parameters {su, si 2 , S13, su, «2i, ■ ■ ■ , s<^}- 

< = n n ( c ™T m w 

u=l m=l 

1, if CW is interested 

w um (23) 

0, otherwise 
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IV. Simulation Results 

In many applications [2]-[4], a 3D cat map is used over a module N as a finite state system shown 
in Eq. (24), where the vector [xt,yt,zt] T and [xt+i, yt+i, zt+i] T denote the discrete spatial coordinates 
at the time t and t + 1, respectively. 



x t +i 
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The period length of this 3D cat map can be defined in Eq. (26), where I is the 3 x 3 identity matrix. 

P N (C 3D ) = arg min { (C 3D )* mod N = I J (25) 



TABLE III: Averaged Period Lengths of 3D Cat Maps 











3D Cat Maps 
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ci D 


C3D 


C3D 


C 3 P D 




r*3D 


C3D 


10 


30.9 


53.7 


8.9 


20.7 


58.1 


59.2 


67.2 


20 


50.7 


85.0 


14.4 


32.3 


88.5 


90.2 


99.2 


30 


142.3 


320.9 


19.8 


55.4 


357.6 


362.2 


442.6 


40 


81.3 


143.6 


24.4 


51.3 


156.5 


154.9 


176.7 


50 


156.7 


274.2 


49.3 


109.3 


292.9 


293.0 


330.4 


60 


217.4 


475.5 


30.0 


73.2 


492.1 


512.3 


589.0 


70 


447.8 


961.5 


46.1 


244.5 


1045.7 


1079.4 


1228.3 


80 


155.3 


276.1 


47.3 


92.2 


300.1 


301.2 


330.4 


90 


315.2 


673.0 


46.8 


111.8 


764.8 


762.1 


929.5 


100 


247.3 


420.3 


77.0 


164.9 


439.3 


439.4 


505.7 



Since any orbit in a finite state system is periodic, the system randomness can be largely reflected 
by averaged period length. And it is desired to have a 3D cat map with longer period lengths [15]. In 
regarding to the randomness of 3D cat maps, we perform the following comparisons on averaged period 
lengths using computer simulations. Specifically, we test six 3D cat maps C| D , C^P, Cff, Cp D , , 
C^? 2 and C\y , and measure the period length of each cat map with a random set of parameters. 
Repeat this experiment 10,000 times and calculate the averaged period length for a 3D cat map under 
the module N denoted as P N (C 3D ) defined in Eq. (26). 

10000 

P N (C^) = Pn(C 3 j D ) (26) 

3=1 
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where C? D is the jth randomly generated C 3D matrix. Simulation results of these averaged period lengths 
for N = {10, 20, ■ ■ ■ , 100} are given in Table III. These results clearly indicate that the new proposed 
3D cat maps C^ i , C^? 2 , and C^? m have much longer averaged period lengths than cat maps [l]-[4]. 

V. Conclusion 

In this letter, we have proposed a new family of 3D cat maps with eight parameters, two parameters 
controlling the cat map spatial configuration and the other six controlling the cat matrix elements. It 
incorporates the conventional 3D cat maps proposed by Lian et al. [1], Chen et al. [2], Liu et al. [3] 
and Pan et al. [4] as special cases. It also outperforms these maps by providing more independent 
parameters and longer averaged period lengths. Both improvements on 3D cat maps are very meaningful 
for enhancing the security of chaotic cryptosystems and image encryption algorithms [l]-[4], [7], [9]- 
[11]. A 3D cat map based system using new proposed 3D cat maps will have a larger key space to 
resist brute-force attacks and a longer averaged period length to resist statistical attacks. The presented 
framework is "universal" and allows extending to higher dimensional cat maps [8], [16]. 
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